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A new approach to dissipative quantum systems modelled by a system plus environment Hamiltonian is 
presented. Using a continuous sequence of infinitesimal unitary transformations the small quantum system 
is decoupled from its thermodynamically large environment. Dissipation enters through the observation 
that system observables generically "decay" completely into a different structure when the Hamiltonian is 
transformed into diagonal form. The method is particularly suited for studying low-temperature proper- 
ties. This is demonstrated explicitly for the super-Ohmic spin-boson model. 
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Dissipative behaviour in quantum physics can be mod- 
elled in a system plus bath framework M . The quantum 
system that one is interested in is described by a Hamilto- 
nian H$ operating on the system Hilbert space 7is, the 
thermodynamically large environment by some Hamil- 
tonian Hb operating on TLb- The full Hamiltonian 
H : TLs^Hb — ► 7~Ls®'Hb of system plus bath is obtained 
by coupling TLs an d 7~Lb with some interaction H$b 



H = Hs <£> lj 
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Most theoretical work starts off by tracing out the bath 
degrees of freedom and then using suitable approximation 
schemes for the time evolution of the reduced density 
matrix of the small quantum system (for a review see 
e.g. Ref. 0). In this letter we present an alternative 
approach that aims at decoupling system and bath with 
a unitary transformation U 
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Here Hs and Hb are modified system and bath Hamil- 
tonians. By carrying out this programme in the manner 
described below this approach is particularly suited for 
studying low-temperature properties of dissipative quan- 
tums systems, thereby being complementary to most 
other approximation schemes. As a specific example we 
demonstrate these ideas for the spin-boson model 
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describing a two-level system coupled to a bath modelled 
by harmonic oscillators. The standard approach to this 
problem is the "Non-Interacting Blip Approximation" 
(NIBA) for the effective action obtained after integrating 
out the bath degrees of freedom |^] . In our approach we 
find low-temperature equilibrium correlation functions 
of the tunneling particle that combine NIBA-results at 
intermediate time scales with the correct long-time be- 
haviour where the simple NIBA fails M. The universal 



Wilson ratio for a super-Ohmic bath put forward in Ref. 
is also obtained. 

Extension of our scheme to other dissipative quantum 
systems is straightforward under the basic assumption 
that Hs has a non-degenerate ground state separated 
by a finite gap from its excited states. Some technical 
details of our method can be found in Ref. ||. 

Two obvious questions arise with respect to the pro- 
gramme (0)-(||). Where is dissipation in Eq. (||) as ex- 
change of energy between system and bath is no longer 
possible? And how can one find a unitary transforma- 
tion U that fulfills the required task? 

We first discuss the second point. The decoupling of 
system and bath is achieved by a method of infinitesimal 
unitary transformations ( "flow equations" ) introduced by 
Wegner in Ref. |(|. A suitable antihermitean generator 
j](£) — —r){ly is chosen and the initial-value problem 



dH{£) 
d£ 
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solved. The parameter £ has dimension (Energy) -2 . 
Eq. (^]) generates a one-parameter family of unitarily 
equivalent Hamiltonians H(£). In the limit £ — > oo one 
attempts to obtain a Hamiltonian H(£ — oo) of the sim- 
ple structure (0). The choice of rj(£) is inspired by renor- 
malization theory: Initially for small £ matrix elements 
corresponding to large energy differences between system 
and bath are decoupled, for large £ one deals with the 
nearly resonant modes. The fundamental problem of (^) 
is that higher and higher order interactions are succes- 
sively generated. So a further condition for rj{£) is that 
the number of additional terms should be small. These 
conditions essentially fix a unique structure of rj{£) 
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All parameters in (g) depend on £ and normal-ordering 
with respect to the noninteracting ground state has been 
introduced. For the construction of generators r\ in a 
general setting see Refs. The only new interaction 

terms generated in the first application of (||) is : a x (b k ± 
b\){b k , ±b\,) :. 

Due to the flow the spectral function J(oj) = 
J2k — w fe) describing the coupling of sys- 

tem and bath becomes ^-dependent too, J(u>, £) — 
J2 k AfcW 2 <5(o, - ui k ), J(uj,£ = 0) = J(lu). We end up 
with the following set of differential equations for the 
couplings in the Hamiltonian by comparing the lhs and 
rhs of M 
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and a differential equation for an uninteresting additive 
term in (j3|). The differential equation for the higher 
normal-ordered interaction term : a x (b k ± b k )(b k / ± : 
is subsequently negelected. This approximation can be 
systematically improved by taking such higher-order in- 
teractions that are also of higher order in the small pa- 
rameters Xk into account one after the other in the hierar- 
chy of differential equations. However, already the above 
approximation leads to a very satisfactory description as 
will be shown. 

As required large energy differences are first decoupled 
in (H|) and small energy differences later. For £ — > oo the 
coupling J(uj,£) vanishes for all o,, in general exponen- 
tially and for u> = A(t — oo) algebraically. Within the 
above approximations one ends up with a system Hamil- 
tonian Hs — — \/2A r a x that is decoupled from the en- 
vironment. Here A r = A(£ = oo) is the renormalized 
tunneling matrix element. By either numerical solution 
of the differential equations or analytical approximations 
one finds A r = cA exp ( — {lo c /K ) s_1 /(2s — 2)) for a 
super-Ohmic bath J(uj) = K 1 ~ s uj s Q(<jJ c — oj) with s > 1. 
Here Ao = A(£ = 0), c is a constant of order 1, ui c some 
high-energy cutoff and K the coupling constant. A r de- 
fines the low-energy scale of the problem in agreement 
with other methods (see e.g. Ref. ||). 

For small £ <C A~ 2 the above procedure is equivalent to 
Anderson's "poor man's" scaling [JtJ with a smooth cut- 
off, except that the high-energy states are not removed 
by integrating them out but by decoupling them using 
unitary transformations (see Fig. 1). By explicitly find- 
ing this unitary transformation no information about the 
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FIG. 1. Schematic behaviour of the effective spectral func- 
tion J(uj,1) for various regimes of the flow equations. 



removed states is lost but contained in the unitary trans- 
formation itself [|| . Scaling approaches have to stop when 
the effective band edge becomes of order the low-energy 
scale of the problem due to divergencies in the renor- 
malization group equations. The flow equations can be 
integrated further since with our choice of /(o,, £) decou- 
pling is with respect to energy differences and does not 
only take place at the effective band edge. 

In the region I > A~ 2 that remains unreachable in the 
"poor man's" approach new features appear that are hid- 
den for smaller £. First of all the flow of the parameters in 
the Hamiltonian becomes negligible and has the universal 
algebraic behaviour £~ 1 / 2 . In contrast the transforma- 
tion of the system observables turns out to be essential 
in this regime. In the crossover region £ w A~ 2 both flow 
of parameters and system observables are important. 

As a specific example for this scenario we discuss the 
symmetrized equilibrium correlation function describing 
the tunneling particle C(t) = \ < {cr z (t), c z (0)} > • In 
order to use the trivial time evolution with respect to 
the Hamiltonian H{£ — oo) we have to transform the 
observable a z as well 



do z (l) 
d£ 
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These differential equations cannot be solved in closed 
form and we have to make an ansatz for the transforma- 
tion of a z {£) 



a z (£) = h(l)a z + a x J2 Xk(£)(b k + b{) 
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where higher normal-ordered terms are neglected. One 
obtains the following differential equations 
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One can prove /i(£ = oo) = if A r lies in the support 
of J(w). Therefore a z (£) in ( |Io| ) decays completely un- 
der the sequence of unitary transformations, which is es- 
sential to see dissipative behaviour with a Hamiltonian 
like (||). In general one can show for a system observ- 
able O ® Is that does not commute with the algebra 
spanned by [Hs <8> 1Lb,Hsb]' If some excitation energy 
from the ground state of Hs lies in the support of the 
spectral function then the observable decays completely 
when transforming from (|l|) to (^]) in the sense that no 
such term of structure O <8> 1b survives f§. This result 
also holds in the zero-temperature limit. 

The one-sided Fourier transform C(u) of the corre- 
lation function C(t) = J °° du> C(uj) cos(wi) can be ex- 
pressed as C(lu) = J^k Xfc(°°) coth(/3wfe/2) 5(ui — uj k ) and 
the Xfe(oo) have to be found numerically. This is simpli- 
fied by the following observation. For £ — ► oo the remain- 
ing spectral function J(u>,£) is strongly peaked around 
A r (Fig. 1). In this limit the exactly solved dissipative 
harmonic oscillator H 



H = Atfb- 
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becomes equivalent to a two-state system as the mean 
occupation number < b > of the dissipative harmonic 
oscillator at zero temperature goes to zero with the width 
of the spectral function. The corresponding set of differ- 
ential equations ([ryq jlljjl^ ) for the dissipative harmonic 
oscillator can be solved in closed form ||. This solu- 
tion is exact as no higher-order terms in the Hamilto- 
nian or the transformation of observables appear. For- 
mally this can be seen by introducing functions Sq(z, £) = 
Ysk^t/iz- uil), Si{z,£) = E fcV ^AxfcA fc /( 

s 2 (zJ) = Y, k xi/(z- 



{Jfy and by showing 
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S 2 (z, oo) = s 2 (z,e ) 
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For the dissipative harmonic oscillator (Q) holds exactly 
without the term 0(£q 1 ). Numerically the correlation 
function C(u>) = — 2lo/tt^s S^w 2 — i0+,oo) has been 
obtained by integrating the flow equations up to some 
£q = (2AA r )~ 2 , A <C 1, thereby obtaining S 2 (z,£ ) and 
then adding the second term from ( |l4] ) that describes a 
dissipative harmonic oscillator with the spectral function 
J(oj,£q). Therefore the resolution of the peak in C(uj) 

— 1/2 

for lj = A r is not limited by £ . By using the analogy 
with the dissipative harmonic oscillator we are restricted 
to low temperatures T <C A r . 
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FIG. 2. Exemplary correlation function for a super-Ohmic 
bath (s = 2, K = 40, u> c = 80). The different curves are ob- 
tained by using Eq. (jL4|) for various values of £o = (2AA r )~ 2 . 
The scale is set by Ao = 1. For comparison the NIBA-curve 
is also shown. The maximum of the NIBA-curve is slightly 
shifted with respect to the flow equations and has been suit- 
ably rescaled in order to identify the peaks. 



Some zero-temperature correlation functions obtained 
in this manner are shown in Fig. 2. The final results vary 
very little with A as long as A < 0.5. This gives a posteri- 
ori justification of our approximations. For £ <§; A,T 2 the 
neglected terms are irrelevant in the usual scaling sense, 
for £ » A~ 2 our equations are closed due to the anal- 
ogy with the dissipative harmonic oscillator. As the final 
results for super-Ohmic baths do not depend on where 
these two parts are matched in the crossover region, it 
is reasonable to argue that the approximations are also 
good for £ w A" 2 |(|. 

The correlation functions obtained for such parameters 
are in good agreement with the NIBA for intermediate 
time scales (Fig. 2). That is for u> « A r the curves are 
well-described by a Lorentzian with a peak at u> = A r 
and half width tt/2 J(A r ) Q. For long times the NIBA 
fails as it predicts an exponential decay at T — [0, 
whereas the long-time behaviour is known to be deter- 
mined by the low-frequency behaviour of the spectral 
function. This is made explicit in the Shiba-relation |l^] 
generalized to super-Ohmic baths in Ref. B (note our 
normalization J °° C(u})dui = 1) 
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CM 
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with the static susceptibility \o- Xo can be extracted 
with a Kramer's-Kronig relation and a fluctuation- 
dissipation theorem xo = 1/2 L C '(w) / 'to dus . The nu- 
merical solution of the flow equations is in excellent 
agreement with the Shiba-relation (Table 1). One ob- 
serves only small deviations that disappear in the limit 
lo c ^> Aq, A r = const, (or uj c /K = const.). 
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TABLE I. Representative results from the numerical 
solution of the flow equations for super-Ohmic baths 
J(w) = K^ s u a @{u c -uj). The scale is set by A = 1. _R< thcor) 
is the Wilson ratio from ([[9]) . Numerical errors for the Shiba- 
and the Wilson ratio are estimated as 2%. 
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to the spin-boson model. The method is an exten- 
sion of "poor man's" scaling j?J as it allows to decouple 
modes below the low-energy scale of the model. Instead 
of renormalization group equations with respect to the 
effective band width we have differential equations with 
respect to energy differences parametrized by I. An es- 
sential new feature as compared to the scaling approach 
is the transformation of the observables once the decou- 
pling has reached the low-energy scale. For dissipative 
quantum systems our method resulted in a complete de- 
coupling of the quantum system from its environment. 
Thereby we have successfully matched formal solutions 
Jl yielding the Shiba-relation and universal Wilson ra- 
tios with the well-established NIBA at intermediate en- 
ergies in one consistent scheme. 



Another interesting low-energy property is the impu- 
rity contribution to the specific heat Cj (T) . This quantity 
is trivial to obtain from a Hamiltonian like (||) . But it fol- 
lows from the difference of two extensive quantities and 
it becomes necessary to discuss the flow equations for the 
bath energies too. That could be ignored before since the 
couplings A/c scale with 1 / yN where N is the number of 
bath modes. One finds 
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and the impurity contribution to the internal energy for 
T <C A r is simply 
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For super-Ohmic baths J(to) oc oj s for small uj, s > 1, 
one derives the following expression for the impurity con- 
tribution to the specific heat Cj 



dEj 
dT 



Cl = sT s T(s + 2)({s + l) 



d£A(£)g(£) (18) 



with g(£) defined from J(u,£) = g(£)io s + O(oj s+1 ). In 
particular the impurity contribution to the specific heat 
is not Schottky-like but scales as T s . A sensitive test is 
provided by the Wilson ratio R s = liniT^o c i{T)/(xo T s ) 
also generalized to super-Ohmic baths in Ref. H 



[1] A. O. Caldeira and A. J. Leggett, Physica 121A, 587 

(1983) ; 130A, 374(E) (1985); Ann. Phys. (NY) 149, 374 

(1984) ; 153, 4 45(E) (1984). 

[2] U. Weiss, Quantum Dissipative Systems (World Scien- 
tific, Singapore, 1993). 

[3] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. 
Fisher, A. Garg and W. Zwerger, Rev. Mod. Phys. 59, 1 
(1987); Erratum Rev. Mod. Phys. 67, 725 (1995). 

[4] M. Sassetti and U. Weiss, Phys. Rev. Lett. 65, 2262 
(1990). 

[5] S. K. Kehrein and A. Mielke, to appear in Ann. Phys. 

(Leipzig) (1997). 
[6] F. Wegner, Ann. Phys. (Leipzig) 3, 77 (1994). For recent 

similar ideas in high-energy physics see also S. D. Glazek 

and K. G. Wilson, Phys. Rev. D 49, 4214 (1994). 
[7] P. W. Anderson, J. Phys. C 3, 2346 (1970). 
[8] In principle it is possible to reconstruct the full unitary 

transformation C/ in (Q) by evaluating ^-ordered products 

of rj(£), but the differential approach is more appropriate. 
[9] F. Haake and R. Reibold, Phys. Rev. A 32, 2462 (1985). 
[10] For an Ohmic bath J(u>) oc uj with intermediate and 

strong coupling the final result for C(u>) depends on A. 

This is due to our choice of rjk, q that generates an IR- 

problem for Ohmic baths || . Hence we restrict ourselves 

to super-Ohmic baths. 
[11] This can also be shown analytically by an approximate 

solution of the set of differential equations, see Ref. . 
[12] H. Shiba, Prog. Theor. Phys. 54, 967 (1975). 



R s = sT(s + 2)C(a ■ 



1) km 

x^0 X s 



(19) 



Wilson ratios obtained from the numerical solution of the 
flow equations can be found in Table 1. Agreement with 
( |l9| ) is excellent in the limit u> c ^> Ao, A r = const. 

Summing up, we have applied a new approximation 
method based on infinitesimal unitary transformations 



4 



